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Synchronization in Networks of Linearly Coupled 
Dynamical Systems via Event-triggered Diffusions 

Wenlian Lu, Member, IEEE , Yujuan Han, Tianping Chen, Senior Member, IEEE 


Abstract —In this paper, we utilize event-triggered coupling 
configuration to realize synchronization of linearly coupled dy¬ 
namical systems. Here, the diffusion couplings are set up from 
the latest observations of the nodes of its neighborhood and the 
next observation time is triggered by the proposed criteria based 
on the local neighborhood information as well. Two scenarios are 
considered: continuous monitoring, that each node can observe 
its neighborhood’s instantaneous states, and discrete monitoring, 
that each node can only obtain its neighborhood’s states at 
the same time point when the coupling term is triggered. In 
both cases, we prove that if the system with persistent coupling 
can synchronize, then these event-trigger coupling strategies can 
synchronize the system, too. 

Index Terms —Linearly coupled dynamical systems, Synchro¬ 
nization, event-triggered diffusions, Continuous and discrete 
monitoring. 

I. Introduction 

YNCHRONIZATION of coupled dynamical systems have 
been widely studied over the past decades CQ-0, which 
can be characterized by that all oscillators approach to a 
uniform dynamical behavior and generally assured by the cou¬ 
plings among nodes and/or external distributed and cooperative 
control. 

In most existing works on linearly coupled dynamical 
systems, each node needs to gather its own state and neighbors 
states and update them spontaneously or in a fixed sampling 
rate, which may cost much. In order to reduce the sampling 
rate of the coupling between nodes, specific discretization is 
necessary. As pointed out by ifTOl . event-based sampling was 
proved to possess better performance than sampling periodi¬ 
cally in time. Hence, iTTl - [|T2l suggested that the event-based 
control algorithms to reduce communication and computation 
load in networked coupled systems. And mm showed 
that the event-based control maintains the control performance. 
Event-based control seems to be suitable for coupled dynami¬ 
cal systems with limited resources and many works addressed 
the event-triggered algorithms. fl3| considered centralized 
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formulation, distributed formulation event-driven strategies for 
multi-agent systems and proposed a self-triggered setup, by 
which continuous measuring of the neighbor states can be 
avoided. I201-I2P studied the stochastic event-driven strate¬ 
gies. flU introduced event-based control strategies for both 
networks of single-integrators with time-delay in communi¬ 
cation links and networks of double-integrators. By using 
scattering transformation, COD investigated the output syn¬ 
chronization problem of multi-agent systems with event-driven 
communication in the presence of constant communication 
delays. In some cases, event-driven strategies for multi-agent 
systems can be regarded as linearization and discretization 
process. For example, as mentioned in the following 

algorithm 

m 

X l (t + 1) = f(x'(t)) + f{x J (t))) (1) 

3=1 

can be a variant of the event triggering (distributed, self 
triggered) model for consensus problem. In centralized control, 
the bound for (t\ +1 —t\) = (tk+i ~tk ) 1° reach synchroniza¬ 
tion was given in the paper fl5l when the coupling graph is 
indirected and in ca for the directed coupling graph. 

Motivated by these works, we apply the idea of event- 
triggered sampling scheme to the coupling configurations to 
realize synchronization of linearly coupled dynamical systems. 
Here, for each node, the coupling term is set up from the 
information of its local neighborhood at the last event time 
and the event is triggered by some criteria derived from 
the information of its local neighborhood. That is, once the 
triggering rule of node is satisfied, the coupling term of this 
node is updated. Thus, the coupling terms are piece-wise 
constant between two neighboring event times. We consider 
two scenarios: continuous monitoring and discrete monitoring. 
Continuous monitoring means that each node can observe its 
neighborhood’s instantaneous information, but discrete moni¬ 
toring means that each node can only obtain the its neighbor¬ 
hood’s information at this event triggered time. As a payoff 
for small cost of discrete monitoring, the triggering events 
happen more frequently than continuous-time monitoring. For 
each scenario, it is shown that the proposed event-triggered 
strategies guarantee the performances of the nominal systems. 

This paper is organized as follows. In sec. II, we propose 
event-trigger coupling strategies to guarantee synchronization 
by employing continuous monitoring. In sec. Ill, we consider 
discrete monitoring. Simulations are given in sec. IV to verify 
the theoretical results. We conclude this paper in sec. V. 
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II. Continuous-time Monitoring 


along © is 


We consider the following network of coupled dynamical 
systems with piece-wise constant linear couplings: 

-jj- = f{x\t)) - cJ^LijTixi ( 4 ) - ^( 4 ))> 

3 = 1 

t\<t< 4 + i, i - !,-•* ,ra. ( 2 ) 

Here, P(£) G M n denotes the state vector of node i, the 
continuous map /(•) : R m —»• R m denotes the identical node 
dynamics if there is no coupling. L = [Lij]™j =1 G R m,m is 
the Laplacian matrix of the underlying bi-graph Q = {V, E}, 
with the node set V and link set E: for each pair of nodes 
i 7 ^ j, Lij = — 1 if i is linked to j otherwise Lij = 0, 
and La = — the graph that we consider in this 

paper is undirected and connected, so L is irreducible and 
symmetric, c is the uniform coupling strength at all nodes, 
and T G R n,n denotes the inner configuration matrix. Let 
0 = Ai(L) < A 2 (L) < • • • < A m(L) be the eigenvalues of L 
with counting the multiplicities. 

The increasing triggering event time sequence {t\}^ =1 (to 
be defined) are node-wise for i = 1, * • • , m. At time t , each 
node i collects its neighbor’s state with respect to an identical 
time point 4 -m with ki(t) = argmaxk'{t k ' < t}. 

For the node dynamics map /, we suppose it belong to some 
map class Quad(P , cP, (3) for some positive definite matrix 
P g R n,n , constant aGR and positive constant /3 > 0, i.e., 


(u — v) T P 


f(u) - f(v) - aP (u - v) 


< —/3(u — v) T P(u — v ) 


( 3 ) 


holds for all u,v G M n . In fact, we do not need this Quad 
condition © holds for all u, v G R n but for a region D c R n , 
which is assumed to contain a global attractors of the coupling 
systems ©. 

We highlight the basic idea behind the setup of the coupling 
term above. Instead of using the spontaneous state from the 
neighborhood to realize synchronize, an economic alternative 
for the node i is to use the neighbor’s constant states at the 
nearest time point t\ until some pre-defined event is triggered 
at time 4 +P then the incoming neighbor’s information is 
updated by the states at t l k+1 until the next event is triggered, 
and so on. The event is defined based on the neighbor’s and 
its own states with some prescribed rule. This process goes on 
through all nodes in a parallel fashion. 

To depict the event that triggers the next coupling time point, 
we introduce the following candidate Lyapunov function: 

1 m m 

V(t) = -xHt)) T P(x\t) 

i= 1 3=1 

= x T (t)(L 0 P)x(t) (4) 

with x = [; x lT (t ),••• ,x rnT (t)] T G R nm and 0 represents 
the Kronecker product. For a compact expression, we denote 
F(x) = [/(x 1 ) T , • • ■ , f(x rn ) T ] T . Then, the derivative of V(t) 


-V(t) J 2 J =2x 1 (t)(L ® P)<j F(x(t)) - a(I m ® T)x(t) 


- [(cL - al m ) ® I ]•!'(/) + ce(t) 


(5) 


where e(t) = [ej ( t ), • • • , e* ( 4 ] T and 

e«(4 = E L v T (*44 - ~ ®*(4 + *44^))) • 
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By assuming / G Quad(P , cP, /?), we have 




l 2 l 


< - 2 fix 1 {L 0 P)x - 2 x t [L(cL - al m ) 0 (PT)]x 
+ 2 cx T (L 0 P)e 

< - 2 fi'x T (L 0 P)x - 2(£ - fi')x T {L 0 P)x 

- 2(cA 2 (L) - a)x T [L 0 (PT)]x + 2 cx T (L 0 P)e (6) 


with any 0 < fi' < fi. Pick a constant v > 0, then 

2x T (L 0 P)e < ux t (L 2 0 P 2 )x + —e T e 

v 

< v\ rn (L)X rn (P)x T (L 0 P)x + ie T e. (7) 

v 

Substitute inequality Q into ©, we have 


jVit) m <-2(J3-F)V(t)-2(cX 2 (L)-a)V(t) 


Denote 


+ [ — + cuA m (L)A m (P)]V r (t) + “ eTe - (8) 


- * < (4) Tp ( a;J '(4 - *44)- 
V J#* 


Then, we have F(t) = Y1T=i 

Denot e || • || the Eucli dean norm, i.e., for any vector £ G M n , 
||£|| = V£? + • • • + £ 2 - For a matrix A G R nXn , the spectral 
norm of A is induced from Euclidean norm, i.e., ||A|| = 
V^max(A T A). Hence, ||A£|| < ||A||||£|| always ho lds, wh ich 
will be used later as default. Moreover, ||x||p = \fx T Px for 
some positive definite matrix P G R n,n . Thus, we have the 
following theorem 


Theorem 1 . Suppose that f G Quad(P , aP, (3) with positive 
matrix P and /3 > 0, cA 2 (L) > a and PT is semi-positive 
definite. Pick fi' < fi. Then either one of the following two 
updating rules can guarantee that system da synchronizes: 

1) Set 4+1 as the time point by the rule 


4+i = max { T > 4 


e i( T )ll < 




cyj A m (L)A m (P) 




(9) 


2) Set 4 +1 os the time point by the rule 


4+i = max {r > t\ : ||ei(r)|| < aexp(-6r)} . ( 10 ) 
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Proof: Noting that c\ 2 (L) > a holds and PT is semi¬ 
positive definite. By (0), we have 

j t V(t) J 2 J < - - P)V(t) + [cu\ m (L)\ m (P) 


-20'}V(t) + -\\e(t)f 

m 

= - 2 (|8 - &)V{t) + J2lcv\m(L)\m(P) 

i=1 
m 

- 2 p']zUt) + J 27 M t ^ 2 (11) 


i=1 


holds for any v > 0 . 
(1). In case of 


< -[2 p - cvX m (L)X m (P)]zf(t) (12) 
c 


for some constant v > 0 , we have 
d 


dt v(t) m <-w-F)v. 

This implies that V(t) converges to 0 exponentially. Note 

max ^[2/3' - cvA m (L)A m (P)] = f— 

^>0 c c A m (L)A m (P) 


(13) 


Then, we take v 
holds. 

(2). In case of 


J 8 ' 


cA m (L) A m (P) 


, which guarantees that © 


\ e iif) || < a exp (— bt ). 


(14) 


we have 
d 


j t V{t) 0 < - 2(/3 - P)V + [cvX m (L)X m (P) 


a 2 cm 


Pick v = 


2/3' 


cXm (L)Xm(P) 


— 2/3']V H-exp (— 2 bt). 


. Then, we have 


dt 


V(t) J 2 J <-2(0-0')V 


arc^mXm(Jj)X 

m (f) 

2v& 


exp (—2 bt), (15) 


which implies that V (t) converges to 0 exponentially. 

■ 

In fact, in Q, if t = t\ but the system does not synchronize, 
then the left-hand term ||e^(f^)|| = 0 and at least one agent 
has positive right-hand term, which means the next inter-event 
interval of one agent must be positive. While in (flOl) . if r = 
t\ but the system does not synchronize, the left-hand term 
equals 0 and the right-hand term is positive for all nodes. 
Hence, the inter-event intervals of all agents are positive. But 
in case the derivative of ei(t) is sufficiently large, the inter¬ 
event interval might tends to 0. Since the dynamics of e*(£) is 
highly related to the property of /(•), towards a lower-bound 
of the inter-event intervals, we suppose that /(•) is Lipschitz 
in the following theorem. 

It should be highlighted that if / is Lipschitz and PT is 
semi-positive definite, then / E Quad(P , aT, /3) with /3 = 


aAi(Pr) 
A m(P) 


x m(P) 
Xl (P) 


. Thus, we have 


Theorem 2. Suppose f E Quad(P , aT, /3) with positive ma¬ 
trix P and [5 > 0 , satisfies Lipschitz condition with Lipschitz 
constant Lf, and there exists some a (possibly negative) such 
that 

(u — v) T P(f(u) — f(v)) > a(u — v) T P(u — v ) (16) 

for all u,v E M n . c\ 2 (L) > a and PT is semi-positive 
definite. For any [3' < [3, any initial condition and any time 
t < 0, we have 

1) With the updating rule at least one agent has 
next inter-event interval, which is lower-bounded by a 
common constant to > 0 . in addition, if there exists 
g > 0 such that zf(t) > gV(t) for all i — 1 , • • • , m and 
t > 0, then the next inter-event interval of every agent 
is strictly positive and is lower-bounded by a common 
constant. 

2) With the updating rule dza the next inter-event interval 
of every agent is strictly positive and is lower-bounded 
by a common constant. 

Proof: (1). Note 

xgt) - x\t) =f{x i ) - f{x i ) - c[(L ® r 

+ c[(L ® T)x(tl)]i, t e [tfc, tfc+i). 

Combining with the facts that / is Lipschitz and V(t) is 
decreasing, we have 

\\x\t)-x\t)\\<L f \W{t)-x\t)\\+c\\T\\V(t\). (17) 

According to 

efit) = y^ LjjT j [x J (s) - x z (s)]ds (18) 

j Jt k 

and inequality (fTTl) . we have 

!Mi)ll <^== E/ll r ll + mc ll r ll 2 ) - 4 ). 

(19) 

And, noting V(t) = YliLi z i(f)’ there exists i* such that 

(20) 

V m 

From the condition (IT6l) . (0 gives 

V > 2 ax T (t)(L (g> P)x(t) - 2 cx T (t)[L 2 ® ( PT) s }x(t ) 

+c2x T (t)(L 0 P)e(t). (21) 

Noting 

u t [L 2 ® ( PT) s ]u < A m (£)^TP« T [L ® (P)]u, 

A\F) 

for all u E R mn , 

2 x T (L ® P)e > —vx T (L 2 <g>P 2 )x——e T e 

> -v\ rn (L)X rn (P)x T (L 0 P)x - —e T e, 

v 

and event gives 

V(t) > wV(t) 
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for all t before the next triggering time, with 

W = 2u-2A m (£) "\ G" -v\ m (L)X m (P) 

A i(+J 

i F 

v cA m (L)A 

m (P)' 

Thus, we have V(t) > V(t l k ) exp (w(t—t l k )). Combined with 
(fl9l) and (l20l) . this implies that for each t < t l k + to, where 
to satisfies 

(L / ||r||+mc||r|| 2 )y^ _ _ ^V^exp(roro) 

TMPj T ° " cVA m (L)A m (P) 

we have the inequality in © holds. Therefore, we have the 
next triggering time should be larger than t l k + to . 

In addition, if zf(t) > qV(t) for all i = 1, • • • ,m and 
£ > 0 , replace (l 20 l) by 

Zi{t) > y/cV[t) (22) 

for allz = 1, • • • , m. Then following the same arguments after 
(l 20 l) . we can conclude that we have the inequality in © holds 
for all t>t\+ to' with some positive to' satisfying: 

(L f \\T\\ +mc||T|| 2 ) y/2^ _ /3'gy/exp(mTo') 

VHP) TO ‘ ~ Cy/\m{L)\ m (Py 

(2). Under updating rule (flOt . By inequality (fl5l) , we get 

V(t) < p exp [—2 min( 6 , ft — /3')t] 

with p = U(0) + Q C +1- Hence, combined with 

GH), this gives 

ll«i(0ll <]jj (i/l|r|| + mc||r|| 2 ) 

x exp [- min( 6 , /? - /3')4](* - 4 )- 

Therefore, ci is guaranteed by the following inequality 


neighborhood but rule rfTOl ) does not need. Therefore, rule 
costs more updating times than m. Moreover, as shown 
by Theorem [2] rule m can guarantee the positivity of the 
intervals to next updating time for all agents but rule 0 can 
only guarantee it for at least one agent at each time or for all 
nodes under some specific additional conditions. 


III. Discrete-time Monitoring 


By the discrete-time monitoring strategy, each node i only 
needs its local neighborhood’s state at time-points k = 

1, 2, • • •. By this way, the design of the next t\ +1 depends 
only on the local states at time t\, other than the triggering 
event ©, which requires the continuous time states. For 
early works, see El, El for reference. 

Consider system ([2]) and the candidate Lyapunov function 
V(x) with its derivative ©. To propose a triggering criterion, 
which depends only on t\ by the criterion © in Theorem Q] 
we need to estimate the bounds of ( x q (t ) — x q (t l k )) — (x l (f) — 
x l (t l k )) for any q,i with L iq 0 and x l (f) — x^(t) for any 
i 7^ f 

First, we estimate the lower-bound of Xi(t) — Xj(t ), which 
satisfies 

d[»’(«)-»>(«)] = [f(xi{t)) _ f{xj{t))] + 0i(4j(t)) _ 0 . ( 4 j(t)) 

by provided the initial values at t k : x l (t k ) and x^(t k ). This 
can be generalized as 

Jw = /( m (4) + 0 «(0) = Mo 23 

If =mt))+0 v(0)=v o . 


Suppose that the solutions satisfy the following inequality: 


(m(<) -v(t)) T P(u(t) 



> Q(t,6,i!>,u 0 ,vo) 


(24) 


2 P 


V Ai(P) 

< a exp (—bt) 


(L/||r|| + mc||r|| 2 ) exp [- min(6, /3 - P')t\](t - t\) 


Since at time t = t k , ei(t ) =0 holds. Based on rule (flOt . the 
next event will not trigger until e*(t) = a exp (—bt). Thus, the 
inter-event intervals t = 4+i - 4 is lower bounded by the 
solution td of the following equation 



(i/liril +mc||r|| 2 ) 


exp [- min(6, /3 - /3')4] t d 


= aexp[-b(r D +t\)\. 


It can be seen that this equation has a positive solution. This 
completes the proof. ■ 


Here q can be regarded as the lower-bound estimation of the 
distance (in P-norm) between two trajectories: 

\\u{t)-v(t)\\p 

= [ [f(u(s))-f(v(s))]ds + (9-i9)t-(u 0 -v 0 ) . 

Jo p 

To specify q, the celebrated Gronwall-Bellman inequality is 
used, which can be verified straightforwardly and described as 
follows: 

Lemma 1. KM 1, B For a nonnegative differentiable function 
x(t) and two continuous functions: a(t ) and (3(t), defined on 
interval [a, b\, if 


x > a(t)x(t) + fi(t) 


Remark 1. The updating rules © and are different but 
closely related to each other in some respects. It can be seen 
from inequalities dZZK and Eli used in the derivation that 
the convergence behavior for 0 might be better than m. 
However, it makes rule © more complicated than f70l) . since 
each agent should receive the message of the states of its 


for all t G [a, b\, then we have 

x(t) > x(a)ex p (/ a(s)di 
+ J a 4(4 exP a a(u)du) ; 


( 25 ) 
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for two continuous functions x(t) and f3(t), and an integrable 
function a(t), defined on interval [a,b\, if /?(•) is nonnegative 
and 


x(t) < a(t) + f f3(s)x(s)ds 
J a 

for all t G [a, b], then we have 

x{t) < a(t) + t a(s)/3(s) exp or P(u)du ^ . (26) 


By this lemma, q is a nonnegative-valued continuous map 
and satisfies (i). p(-, 6 , 6 , ito, wo) = 0 ; (ii). p( 0 , •, •, uo, uo) = 0 . 
For example, assuming that condition (fl6t holds, we have 


J t l{u{t) - v(t)) T P(u(t) - v(t))\ 

=2 (u - v) T P[f(u) - f(v) + 6 -ft] 

>2 a(u — v) T P(u — v) — fi(u — v) T P(u — v ) 

- -{e-$) T p(e-$) 

l 1 


for any p > 0. By the Gronwall-Bell inequality (1251) . we have 


(u(t) - v(t)) T P(u(t) - v(t)) 

> exp [( 2 cr - n)t](u 0 - v 0 ) T P{u 0 - Vo) 

2j(T jjj ^ 


which is positive for a small interval of t, starting from 0 , and 

Uo V 0 . 

It can be seen that holds for a large class of functions 
/(•). For example, if there exists some a G M such that 

{P^(x)} s > aP (27) 

for all x G I n and some a G R, then 


p's generator 


Input 


Integrators 


Output 



Fig. 1. p generator. 


displacements of two trajectories of (l23l) with respect to their 
initial locations, 

||(u(t) -u 0 ) - (v(t)-v 0 )|| 

{f(u(s)) - f(v(s))]ds +{6-tf)t . 

For example, if /(•) is Lipschitz (on the two trajectories): 

|| f(u(s)) - f(v(s))W < Lf\\u(s) - v(s)|| for all s > 0 , then 
we have 

||(u(t) -u 0 ) - (y{t) -u 0 )[| 

<Lf ||(u(s) - u 0 ) - (u(s) - v 0 )\\ds 
Jo 

+ (||(9-^|| +L f \\u 0 -v 0 \\) t. 

By the Gronwall inequality (l26l) . we have 



(u - v) T P[f(u) - f(y)] = f (u- v) T 

Jo 

{P^r~(\(v — u) + v)} s (u — v)d\ > a(u — v) T P(u — v ) 
ox 

which implies (IT 6 T) . In particular, if the Jacobin of /(•) is 
bounded, namely, \\df/dx\\ > Jf for some Jf > 0 , then we 
have 

{p%{v)Y > ~Jf\\P\l 

which implies a = —Jf x^ P \p) G3- 

Second, we consider the differential equations (l23l) and sup¬ 
pose that the solutions of (l23l) satisfy the following inequality: 


\\(u(t)-u 0 ) - (v(t) - Vo)II 
^ (llg-tfll+^lluo-uoll) _ (29) 

L f 

It can be seen that the upper-bound equals to zero if t = 0. 

It can be seen that the estimation of p and q substantially 
depends on the form of /(•). There might not be a unified 
approach to give precise estimation for general /(•) but might 
be done case by case. Therefore, an efficient but cost way 
is to use integrators that simulate the node dynamics of ii = 
f(u) + Q as the generators of p and q. These generators are 
independent of the states of the nodes and so parallel to the 
networked systems. Figs. Q] and [2] show the configurations of 
the generators of p and q respectively. 

Let 


II (u(t) - u 0 ) - (v(t) - u 0 ) || < p(t, 0,i9,u o ,v o ), (28) 

where p is nonnegative-valued continuous map that depends 
on the node dynamics map /(•), the initial value uo,vo and 
inputs and satisfies p(0,= 0. Geometrically, p 
is an upper-bound estimation of the difference between the 


e k = - ^(4)). 

3=1 

m 

1=1 
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q's generator 


Input 


Integrators 


Output 



by the continuous dependence of the parameters in the system 
0 , m has indeed a positive maximum £. 

Secondly, each node i needs to know the states of itself 
and its neighbors. In details, when one node j is triggered, it 
sends off its new coupling terms, 6j (t J k ^), to all its neighbors 
for their updating the estimation of f k for their next updating 
times. 

Remark 3. In case that x J ( t k ) = x l {t k ) for all neighbors j of 
node i, both left and right sides equal zero, which might lead 
to a Zeno behavior. To avoid the Zeno behavior, we provide 
a triggering event, which depends only on t k by the rule m 
in Theorem [7] Here, we only need to estimate the bounds of 
(. x q (t ) - x q (t[)) - (x z (t) - x l (t l k )) for any q , i with L iq 0. 
Note 


Fig. 2. q generator. 


(x q (t) - x q (t[)) - (x l (t) - x l (t[)) 



f(x q (s))-f (Vis))+ 91^-61 


ds. 


Then, for t >t k , we have 

ll e *( r )ll < Z,{~ Liq ^ T \\ p ( T - t k’ 9 k’ 9 k q (T)’ xt ( t l)’ x 9 ( t l))’ 

qy£i 

Zi(r) > jYy-LijWiT-t^ei,el. (T) , V(t% ), Xi (4)). 

V 3& 

With these assumptions and Theorem [U we have the fol¬ 
lowing result. 


Combing with \\6 k — / J| < M, where M > 0 is some 

constant, we suppose that the solutions of (E3 satisfy the 
following inequality: 

II (u(t) - M 0 ) - (v(t) - M Q )|| < p 1 (t,u 0 ,v 0 ). (31) 

Then, for t >t\, we have 

ikiWii < (*{))■ 

j^i 


Theorem 3. Suppose that f e Quad(P, aY, /3) with positive Theorem 4. Suppose that f e Quad(P, aT, /3) with positive 
matrix P and 0 > 0 such that PT is semi-positive definite, matrix P and 0 > 0 such that PT is semi-positive definite. 
c\ 2 (L) > a. For any positive 0' < 0, set by: c\ 2 (L) > a. For any positive 0' < 0, set inter-event interval 

a - -{<: * . „ ( £: E( _ Ilj)l|r|lwtt , I . (4W(4)) 


< 


0 ' 


C\/XZJX)Xn{P) . 

x\ti),x j (ti)) 


J2(-L iq )e\ceiel 


< ae-xp[-b(£ + f k )] \. 


(32) 


(30) jj le even f times {t\} are set by the following scheme: 


The event timing {t\} are set by the following scheme: 

1) Initialization: t l Q = 0 for all i = 1, • • • , m; 

2) For node i, set via its neighbor's and its own current 
states and diffusion by m; 

3) If one of its neighbors, for example, j, triggers at t = 
^k'+i (k' ^ th e latest event at node j before t), then 
replace 6 J k , by 0 k , +1 in <E3 and go to Step 2; 

4) Let t k+1 sat t k + f k , the event triggers at node i by 
changing t k in d 2 ]) to t k+1 . 

Then, system 0 synchronizes. 

The proof of this theorem can be derived by using in 
Theorem Q] directly. 

Remark 2. Firstly, if node i has a neighbor j satisfying 
xi (t k ) 7 ^ t ^ ien %k i n (E3 is well defined. In fact, in 

POl) , if f = 0, the left-hand side of the inequality in POl ) 
equals zero while the right-hand side is nonzero. Therefore, 


1) Initialization: t z 0 = 0 for all i = 1, * • • , m; 

2) For node i, search f k via its neighbor's and its own 
current states by 03; 

3) Triggers node i by changing t k in 0 to t k+1 = t l k -\-f, k >- 
Then, system 0 synchronizes. 

This theorem can be derived by using (flOl) in Theorem Q] 

It should be highlighted that under rule (l32i every node 
does not need to know the coupling terms of neighbors 
anymore and the inter-event intervals have a lower-bound. 

In fact, in oi, if £ = 0 , the left-hand part equals zero 
while the right-hand is nonzero. Therefore, according to the 
continuous dependence of the parameters in the system ( 0 , 
(l32l) has indeed a positive maximum argument ££. 

Similar to Theorem [21 we have 

Theorem 5. Suppose that f E Quad(P , aT, (3) with positive 
matrix P and (3 > 0, satisfies Lipschitz condition with Lips- 
chitz constant Lf, and there exists some a (possibly negative) 
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such that m holds for all u,v G W 1 . cA 2 (L) > a and PT 
is semi-positive definite. For any (3' < ft, any initial condition 
and any time t < 0 , we have 

( 1 ) under the updating rule dm there exists to > 0 such 
that there exists at least one agent k E {1, • • • , m} such 
that the next inter-event interval is strictly positive and 
has the lower-bound to; in addition, if there exists q > 0 
such that z?(t) > qV ( t ) for all i = 1 , • • • , m and t > 0, 
then the next inter-event interval of every agent is strictly 
positive and has a common positive lower-bound. 

(2) suppose /(•) is Lipschitz with constant Lf. Then, under 
the updating rule m, the next inter-event interval of 
every agent is strictly positive and has a common lower- 
bound t' d . 

Remark 4. In comparison with the continuous-time moni¬ 
toring, the discrete-time monitoring works well particularly 
when the states of nodes cannot be monitored spontaneously. 
Generally speaking, the main difference between these two 
monitoring strategies is that continuous-time monitoring de¬ 
termines the next updating time in an on-line way, based 
on the spontaneous information of states of nodes. Instead, 
the discrete-time monitoring predicts the next updating time. 
Therefore, the discrete-time monitoring costs less for collecting 
state information than continuous-time monitoring. However, 
as a trade-off, it needs more calculations in predicting the next 
updating time, as mentioned in dm or dm. 



Fig. 3. Topology of the graph of the coupled system and the pinned node. 



IV. Examples 


In this section, we present two examples to illustrate the 
theoretical results. The system is an array of 10 linearly 
coupled Chua circuits with the node dynamics 


m 


p * (-Z1 + Z2- g^Zx)) 
Zl - Z 2 + z 3 


-q*z 2 


(33) 


where g(zx) = mi* 2:1 + 1 / 2 *(mo — mi)*(|zi+ 1 | — \ zi — 1 |), 
with the parameters p = 9.78, q = 14.97, mo — —1.31 
and mi = —0.75, which implies that the intrinsic node 
dynamics (without diffusion) have a double-scrolling chaotic 
attractor El. The coupling graph topology is shown in Fig. 
El L is picked as the Laplacian of the graph where each 
link has uniform weight 1. Then, the largest and smallest 
nonzero eigenvalues equal to A 2 = 0.8363 and A m = 7.3484 
respectively. Let P = T = / 3 . To estimate the parameter (3 in 
the Quad condition, noting the Jacobin matrices of / is one 
of the following 


Ax 

A 2 


-2.445 9.78 0 

1 -1 1 
0 -14.97 0 

3.0318 9.78 0 

1 -1 1 
0 -14.97 0 


then we can estimate f3' = a — A max ((A 2 ) s ) = a — %, where 
X = 9.1207 is the upper-bound of the largest eigenvalues of 
the symmetry parts of all possible Jacobin matrices of /. 


Fig. 4. Variation of /3'/(cy/\ m ) under continuous-time monitoring for 
synchronization with respect to the coupling strength c. 


The ordinary differential equation © is numerically solved 
by the Euler method with a time step 0.001 (sec) and the time 
duration of the numerical simulations is [ 0 , 2 ](sec). 

A. Continuous-time monitoring 

First, we consider the updating rule © . According to cA 2 > 
a, where A 2 (L) is the smallest eigenvalue of L, except the 
unique zero eigenvalue, the supremum of the term f3'/ {cy/Xff) 
is estimated as follows: 

13' cA 2 (L)-A max (A*) ^ A 2 (L) 

SU.J/) / 1 1 / 

C V (L) fy A m (£) VA m (L) 

as c -+ 00, by picking a = cA 2 (L). Fig. |4| shows the variation 
of f3’! (c^/A m (L)) with respect to c. In this example, we pick 
c = 20.3281, which implies /3'/ (c^/A m (L)) = 0.1450. We 
employ the updating rule © in Theorem □ Rg. a presents 
the dynamics of each component of the 10 nodes and show 
that the coupled system © reaches synchronization. Fig. |9] 
shows that V ( t ) decreases with respect to time and converges 
toward zero as time goes to infinity. 

Second, we employ the updating rule (flOl) . We take the same 
value of c as above and a = 0.5, b = 0.5. Fig. | 6 ] presents the 
dynamics of each component of the 10 nodes and show that 
the coupled system © reaches synchronization. Fig. [9] shows 
that V ( t ) decreases with time and converges to zero as time 
goes to infinity. 
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Fig. 5. Dynamics of components of the coupled system <3 under the event 
triggering rule {9j. 





Time(Sec) 


Fig. 6. Dynamics of components of the coupled system J3 under the event 
triggering rule i fTOl . 





Fig. 7. Dynamics of components of the coupled system $2) under the event 
triggering rule Ool 


B. Discrete-time monitoring 

First, we employ the updating rule (l30t . The term 
F/{cy/Xfh) can be directly derived from the arguments above. 
We pick the same c and then F / (cy / Kf) is the same as above. 
We employ the event trigger algorithm (l30l) in Theorem [3] Fig. 
|7] presents the dynamics of each component of the 10 nodes 
and shows that the coupled system (0 reaches synchronisation. 
Fig. 0 shows that V(t) decreases with time and converges to 
zero. 

Second, we consider the updating rule d32b . We pick the 
same c as above and a = 0.5, b = 0.5. Fig. [8] presents the 
dynamics of each component of the 10 nodes and show that 
the coupled system 0 reach synchronization. Fig. 0 shows 
that V (t) decreases with respect to time and converges toward 
zero as time goes to infinity. 

C. Performance comparison 

In comparison, we consider the original linear coupled 
system as follows: 

rlrUt) JIL 

—= /(z’W) - (t) - x\t )), (34) 

i =1 

for i = 1, • • • , m. By the same setups of model and numerical 
approach as above, its performance in terms of converge rates 
of V ( t ) is shown similar with those of event-triggered rules 
© and as comparatively shown by Fig. 0 As for the 
performance of rules m and (ED, since the exponential 
convergence rates are pre-designed, as shown by (ITOt and (l32l) . 
it is not surprising that their convergence rates are not as good 
as However, their updating times of these rules are much 
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Fig. 8. Dynamics of components of the coupled system {2} under the event 
triggering rule {32). 



Fig. 9. Dynamics of V(t) for systems with continuous updating, continuous 
monitoring under rules (f9t.(fl0). discrete monitoring under rules {30} , {32). 

less than those of rules @ and (1301) . as comparatively shown 
in Figs. [T0(a)lfT0(b) 

V. CONCLUSION 

In this paper, we employed event-triggered coupling con¬ 
figurations to realize synchronization for linearly coupled dy¬ 
namical systems. We studied both continuous monitoring and 
discrete monitoring schemes: continuous monitoring scheme 
means that each node collects its neighborhood’s instantaneous 
state, and discrete monitoring scheme means that each node 
obtains its neighborhood’s states at the event triggered time. 
The event-triggered rules were proved to perform well and 



(a) under updating rules {9} and ( 1301 



Nodes 

(b) under updating rules {lOl and {32) 

Fig. 10. Histogram of triggering times of each node in [1, 2](sec). 

can exclude Zeno behaviors, as proved for some cases and 
illustrated by simulations. We showed that there are trade-offs 
between better performance in terms of fast convergence and 
less updating time slots, and between more cost in observation 
of states and more computation load of predicting next updat¬ 
ing times. One step further, there are a few issues, including 
how to estimate the number of updating time slots and its 
dependence on the parameters in the rule and the structure of 
network structure, which merits the future research. 
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